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Abstract 

In this paper, we analyze the eigenfunctions of the edge-based Lapla- 

Qcian on a graph and the relationship of these functions to random walks on 
. the graph. We commence by discussing the set of eigenfunctions supported 

t/5 at the vertices, and demonstrate the relationship of these eigenfunctions 

1 ^ I to the classical random walk on the graph. Then, from an analysis of 

functions supported only on the interior of edges, we develop a method 
T-H for explicitly calculating the edge-interior eigenfunctions of the edge-based 

^ Laplacian. This reveals a connection between the edge-based Laplacian 

and the adjacency matrix of backtrackless random walk on the graph. The 
edge-based eigenfunctions therefore correspond to some eigenfunctions of 
\l the normalised Hashimoto matrix. 
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> 1 Introduction 



The traditional discrete graph Laplacian operator .2 has proved to be a useful 
tool in the analysis of graphs and has found application in a number of areas. 
For example, the heat kernel, which is derived from the graph Laplacian, has 
been used to define graph kernels piT[|12j in the machine learning literature. Sun 
et al[Tn] used the heat kernel on the mesh representing a 3D shape to create a 
heat kernel signature for describing shape. In [T], Aubry et al used wave-like 
solutions of Schrodinger's equation to construct an alternative shape descriptor, 
referred to as the wave kernel signature. The graph Laplacian was used by 
Coifman and Lafon[3] for dimensionality reduction of data. There are many 
other applications graph Laplacian in the literature. 

* Corresponding author: Richard.Wilson@york.ac.uk 

t Edwin Hancock is supported by the Royal Society under a Wolfson Research Merit Award. 



1 



In [31 [B] , Friedman and Tillich developed a calculus on graphs which pro- 
vides strong connections between graph theory and analysis. This approach has 
a number of advantages. It allows the application of many results from anal- 
ysis directly to the graph domain, and opens up the use of many new partial 
differential equations on graphs. As an example, they define a wave equation [7] 
which has a finite speed of propagation, in contrast to the usual wave equation 
on a graph. 

In the graph calculus of Friedman and Tillich, the graph is given a geometric 
realization by associating an interval with each edge of the graph. Functions 
may therefore exist both at the vertices and on the interior of edges. From 
this starting point they develop a divergence and, most importantly, a graph 
Laplacian. This type of Laplacian has found application in the physics literature 
where the interpretation is as the limiting case of a "quantum wire" [9l [Ml [11] . 

The graph Laplacian consists of two parts; namely a vertex-based Laplacian 
and an edge-based Laplacian. Friedman and Tillich also demonstrate that for 
edgewise-linear functions the edge-based Laplacian is zero and the graph Lapla- 
cian reduces to the traditional discrete graph Laplacian. On the other hand, for 
functions where the vertex-based Laplacian is zero, they obtain the edge-based 
Laplacian only. This results in a setting which is substantially different from the 
traditional approach. In the remainder of this paper, we will concern ourselves 
with the edge-based Laplacian. 

While Friedman and Tillich find the eigenvalues of the edge-based Lapla- 
cian, and give some of its eigenfunctions explicitly, they do not give a method 
for computing the entire eigensystem of the graph. In this paper, we demon- 
strate the relationship between the simple eigenfunctions and the random walk 
on the graph. We then give a method for explicit calculation of the remain- 
ing eigenfunctions. This analysis reveals a link between the remaining eigen- 
functions and the backtrackless random walk on the graph, which in turn is 
linked to some properties of the quantum walk on a graph[4j and its Ihara zeta 
ftmction[l3]. Because of this link, the edge-based eigenfunctions correspond 
naturally to eigenfunctions of the normalised adjacency matrix of the oriented 
line graph (the normalised Hashimoto matrix^). 

The remainder of this paper is organized as follows. In Section 2, we briefly 
introduce the formalism of Friedman and Tillich. In Section 3, we detail the 
calculation of the vertex-supported eigenfunctions of the edge-based Laplacian 
and describe their relationship to the random walk. Section 4 contains our main 
result; we provide a method for calculating the remaining eigenfunctions, an 
approach which is linked to the backtrackless random walk. Together with the 
eigenfunctions detailed in Section 3, this completely describes the eigensystem 
in terms of finite matrices on the graph. Finally in Section 5 we give some 
conclusions and suggest future directions of research. 
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2 The Edge-based Laplacian 



In this section, we give brief details of the formahsm of Friedman and Tinich[S]. 
Let G = {V,E) be a graph with a geometric reahzation Q. The geometric 
reahzation is the metric space consisting of the vertices V and a closed interval 
of length If. associated with each edge e & E. The graph has a (possibly empty) 
boundary set dG C G. We assume that the boundary is separated, i.e. that each 
boundary vertex is incident on only one edge. The graph G—{ V,E} excludes 
boundary vertices and any incident edges G= G\dG. We associate an edge 
variable with each edge. Let e = be an edge with interval variable Xe- 

The edge variable equals zero where the edge meets vertex u and equals one 
at vertex v. The start and end vertices are determined by assigning an arbitrary 
orientation to each edge. 

Definition 2.1 (Friedman and TillichjS]). A vertex measure, V is a measure 
supported on the vertices with V{v) > for all v € V. For our purposes it 
suffices to take V{v) = 1 for all v Cz V. An edge measure, £ is a measure 
supported on the interior of edges. £{v) =0 for all v V and the restriction to 
the interior of edges is the Lebesgue measure. 

Let / be a function defined on the graph (on both edges and vertices). We 
take f{u) to mean the value of / at vertex u and /(e, Xe) to mean the value of 
/ at position Xe along edge e. Since we have different volume measures on the 
edges and vertices, we must take care in dealing with integrating factors since 
they are different on edges and vertices. We use dV for the vertex integration 
factor and d£ for the edge integration factor. As a result, we have a two-part 
Laplacian: 

A/ = AvfdV + AsfdS, (1) 

where Ay is the vertex-based Laplacian and A^; is the edge-based Laplacian. 
Since graph Laplacians are usually given as positive definite operators, the edge- 
based Laplacian is minus the usual calculus Laplacian 

Ai^/--Vcalc-V/. (2) 
The vertex-based Laplacian turns out to be 

^vf^^^Y.^e,.-yfUv). (3) 

Here Ug^y is the outward-pointing unit normal. In other words, for an edge (a, b) 
it points from a to 6 at the vertex b, and from & to a at the vertex a. 

Definition 2.2. A function is said to be edge-based ii Ay/ = 0. For edge-based 
functions, the Laplacian consists of only the edge-based part. A/ = AEfdE 

For a function / to be edge-based, the following condition applies. 

^(-l)i---^'V/(e,Xe,„) = OVz;, (4) 

or in other words, the sum of the outward-pointing gradients must be zero. 



3 



3 Vert ex- supported Edge-based Eigenfunctions 



For the remainder of this paper, we assume that the edge lengths on the graph 
are equaL Friedman and Tilhch[T demonstrate the connection between the 
eigenfunctions of the edge-based Laplacian and the eigenvectors of the row- 
normahsed adjacency matrix. This follows directly from the observation that 
As is essentially the familiar Laplacian of calculus and therefore admits eigen- 
vectors of the form /(e, Xe) = C{e) cos(wXe -I- B[e)) where lj is the frequency of 
the eigenfunction, corresponding to an eigenvalue of w^. The eigenfunction is 
edge-based and so applying condition Q to the eigenfunction gives 

fju) - f{v) cos uj ^ ^ 
^ sinw ■ ^ ' 

for any {w, /} for A^;, when oj is not a multiple of tt. We call the frequency of 
the eigenfunction and {w, /} an eigenpair while noting that the corresponding 
eigenvalue is w^. 

Definition 3.1. A principal eigenpair is an eigenpair of A^; with < lo < 2Tr 
where w is the smallest magnitude frequency of a sequence of eigenpairs of the 
form {(jj + 2mT, C(e) cos[(w -I- 2nn)xe + -B(e)]}, n € N with the same coefficients 
B{e) and C(e). 

Let the vector g 7^ be the restriction of an eigenfunction / to the vertices, 
taken in some particular order, i.e. = f{u). All eigenfunctions in the same 
sequence as / have a vertex restriction equal to g. 

The row-normalised adjacency matrix A for interior vertices v £V is given 

by 

^ = (6) 
where A is the usual graph adjacency matrix. 

Theorem 3.2 (Friedman and Tillich 7J). Let Q he the geometric realization of 
a graph and let A he its row-normalized adjacency matrix. Each eigenvalue A 
of A., with A ^ 1} corresponds to two principal frequencies of Q 

cos~^ A and 27r — cos^^A, (7) 

each with the same multiplicity as A. The corresponding eigenvector g is the 
vertex restriction of the sequence of eigenfunctions based on the principal eigen- 
value (jj. 

Equation ([5| allows us to determine the eigenfunctions. For a principal 
frequency lu the principal eigenfunction is /(e,Xe) = C(e) cos(_B(e) -f wxg) with 

C{ef = .gg+gg-2g.g.cos^^ 
sin u) 

tanB(e) = d.cosu-^g.^ ^ 
gy sm uj 
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There are two solutions to Equations ^ and ^ which are {C(e), -Bo(e)} or 
{—C{e), Bo{e) + tt} but both give the same eigenf unction. The sign of C{e) 
must be chosen correctly to match the phase, i.e. so that C(e) cos(i3(e)) = gu- 
Since B and C are uniquely determined, this comprises all the eigenvectors of 
this form. 

Friedman and Tillich give a pair of frequencies as above for each eigenvalue 
of the row-normalized adjacency matrix, but we note that the eigenpair with 
uj = 2tt — cos~^ X actually corresponds to the same sequence of eigenpairs as 
those with w = cos^^A, but with negative values of n. There is therefore a 
single sequence for each eigenvalue with 

UJ = cos^^A, (10) 
f{e,Xe) = C{e)cos[B{e) + {uj + 2nn)xe],neZ. (11) 

The value A = 1 is always an eigenvalue of A. This corresponds to a principal 
frequency of co = in A^; and therefore the corresponding eigenfunction is 
f{e,Xe) = C(e) cos(B(e)) which is constant on the vertices. If dG 7^ then 
this eigenfunction must be zero everywhere. If dG ~ then we obtain a single 
principal eigenpair with uj = and /(e, Xg) — C . 

The value A = — 1 will be an eigenvalue of A if G is bipartite. If dG ^ 
then the eigenfunction must be zero on the vertices. We defer the evaluation 
of eigenfunctions not supported on V for the next section. If dG = then a 
single principal eigenpair exists with lj — t: and f{e,Xe) = C cos{TTXe). The 
eigenfunction alternates in sign between the two partitions of the graph. 

This comprises all principal eigenpairs which are supported on the vertices 
(g 7^ 0)j7]. The eigenfunctions supported on the vertices are therefore directly 
determined by the eigensystem of A. 

3.1 Random Walks and Line Graphs 

The line graph LG{G) = {Vi, Ei) of a graph G is a graph construction in which 
we replace the edges of G with the vertices of LG(G) as follows. Firstly we 
create the symmetric digraph SDG(G) of G by replacing each undirected edge 
with a pair of oriented edges. Each oriented edge of the SDG then becomes a 
vertex of LG(G). These vertices are connected if the head of one oriented edge 
meets the tail of another. The reverse pair of oriented edges are connected, i.e. 
((u, w), (f , u)) is an edge in the LG. 

Vi = {{u,v)eE{SDG)}, 

El = {((m, v), (f , w)), (u, v) e E{SDG), {v, w) e E{SDG)}. 

The oriented line graph OLG(G) = {Vo,Eo) is constructed in the same way 
as the LG except that reverse pairs of oriented edges are not connected, i.e. 
((u, v), {v, u)) is not an edge. The vertex and edge sets of OLG(G) are therefore 

Vo = {iu,v)eE{SDG)}, 
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Figure 1: Left: A graph. Middle: Its symmetric digraph. Right: The line graph 
and the oriented line graph; the dotted edges exist in the line graph but not in 
the oriented line graph. 



Eo = {((u, w), {v, w)), {u, v) e E{SDG), (w, w) e E{SDG) -.u^w}. 

Figure [T] illustrates these concepts. A random walk on the vertices of LG 
represents the sequence of edges traversed in a random walk on the original graph 
G. Similarly, a random walk on the OLG represents the sequence of edges in 
a random walk on G where backtracking steps are not allowed (a backtrackless 
walk) . 

Proposition 3.3. Let A be the row-normalised adjacency matrix of G and U 
be the row-normalised adjacency matrix of the line graph of G. Each eigenpair 
{A,g} of A corresponds to an eigenpair {/^,h} o/U with 

M - A, (12) 

Proof. We may write the row-normalized adjacency matrix of the LG as 

A A 

f"T — '^'^ ^^^^"^ _ A A s (^A\ 

We have 

^ ^ Uuv jWx^wx — ^ ^ -^uv-^wx^vw-^wxQxi (■^^) 
w,x w.x 

(16) 

X 

= AAu^gt, = y^huv (17) 

□ 

The vertex-supported eigenfunctions of the edge-based Laplacian are there- 
fore determined by the structure of the random walk on the graph. As we shall 
show later, the remaining eigenfunctions are determined by the structure of the 
backtrackless random walk. 
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4 Edge-interior Eigenfunctions 



We now proceed to our main result. All remaining eigenfunctions of the edge- 
based Laplacian are zero on the vertices of G and therefore must have a principal 
frequency of u G {tt, 27r}. We start with the case lu = tt. 

Proposition 4.1. Principal eigenfunctions of the edge-based Laplacian with 
principal frequency uj — n and which are zero on the vertices of Q are of the 
form /(e, x^) — C{e) cos(-| + vrXe), e GiS with 

C(e) = Vv e V. (18) 

Proof. Since the boundary is separated, the eigenfunctions must be zero on 
any boundary edge and any edge incident on a boundary vertex. As a result, 
we may concern ourselves only with G. The eigenfunction f{e,Xe) is zero at 
both vertices incident on edge e, giving values for B{e) of B{e) £ {7r/2, 37r/2}, 
which both give the same eigenfunction (with a different sign for C(e)). We 
may therefore take B{e) = tt/2. The gradients at either end of the edge are 
V/(e, Xe^O) = -7rC(e) sin(f ) and V/(e, x^ = I) = -7rC(e) sin(^). Applying 
condition ^ to this eigenfunction, we obtain 

J2Cie) = 0. 

□ 

Hence in order to find eigenfunctions of this type, we must find a set of 
coefficients attached to the edges which sum to zero at every vertex. 

Proposition 4.2. The principal eigenfunctions of the edge-based Laplacian with 
principal frequency w = 27r and which are zero on the vertices of Q are of the 
form /(e, Xf.) — C(e) cos(|^ + 2TTXe), e £E with 

^(-l)i-^'=-C(e) = OVw e V. (19) 

Proof. The proof is essentially the same as for the previous proposition. However 
the gradients at either end of the edge are V/(e, = 0) = — 27rC(e) sin(|) and 
V/(e, a^e = 1) = -27rC(e) sin(^). Applying Condition (jij), we obtain 

^(-1)1— C(e)=0. 

□ 

We may interpret this condition as follows: Consider each undirected edge 
of G as a pair of directed edges. Associate a value C(e) with each directed edge 
in the direction of increasing Xe (i.e. from Xe — ^ Xe = 1) and a value — C(e) 
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with the reverse edge. Then the condition above is that the sum of incoming 
directed edges at a vertex must be zero. 

We may write these conditions for both principal frequencies in the following 
form. Let = ±C(e) be the values associated with each edge, as described 
above. The sign is always positive for uj = n and alternates in sign depending 
on the direction of edge traversal for uj = 27r. Then we may form a matrix W 
with elements Wuv = Wuv^uv The conditions for uj = tt are then 

Wl = 0, (20) 

W = W"^, (21) 

and for uj = 2tt they are 

Wl = 0, (22) 

W = -W^, (23) 

where 1 is the vector of all-ones. 

The oriented line graph of G (OLG(G)) represents the structure of a back- 
trackless random walk on G in the sense that a random walk on the vertices of 
OLG(G') generate a sequence of edges which can be traversed in a backtrackless 
random walk on G. In [13] we demonstrated that the adjacency matrix of the 
oriented line graph (T) is equal to the positive support of a quantum walk on 
the graph G. This matrix is also related to the Ihara zcta function of the graph 
since it is equal to the Perron- Frobenius operator on OLG(G') and therefore can 
be used to calculate the Ihara zeta function Zq {u) using Zq^u) = det(I 
We now show that the eigenvectors of T corresponding to eigenvalues of A = ±1 
determine the structure of the edge-interior eigenf unctions. 

Theorem 4.3. Let T be the adjacency matrix of OLG{G) and s be an eigen- 
vector of T with eigenvalue A = 1. Then s„„ = —s^u and '^y^Suv = 0, and s 
provides a solution for W in the case of uj — 2tt. Similarly, if X = —1 then 
Suv = and J2u ^uv = 0, and s provides a solution for W in the case of 
UJ = n. 

Proof. In [13] we demonstrate that s^v — ^uvWuv is an eigenvector of T if 
X]„^«t)W„„ = and either w„v = -Wyu or Wuv = Wvu- If Wuv = ~Wvu the 
eigenvalue is A = 1, and if Wuv = Wvu then the eigenvalue is A — —I. These 
are precisely the solutions for W (and for C(e)). Since the eigenvectors with 
A = ±1 span the space of possible solutions, we obtain \E\ — \V\ + 1 linearly 
independent solutions for A = 1 and — |F| linearly independent solutions for 
A = — 1 which are all the available solutions according to [7] □ 

The structure of the eigenfunctions which are not supported on the vertices 
is therefore determined by the eigenvectors of the backtrackless random walk 
on the graph G. 
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5 Conclusions 



We have analyzed and completely determined the eigenfunctions of the edge- 
based Laplacian and given explicit forms for all of these eigenfunctions. Our 
analysis provides a method of computing the eigenfunctions which are zero on 
the vertices from the eigenvectors of the oriented line graph. We demonstrate 
the connection between the eigenfunctions and both the classical random walk 
and the backtrackless random walk. The eigensystem of edge-based Laplacian 
contains eigenfunctions which are related to both the adjacency matrix of the 
line graph of G and the adjacency matrix of the oriented line graph of G. 

As noted by Friedman and Tillich[7j, this approach is closer to traditional 
analysis than the usual discrete graph Laplacian. In particular it allows us to 
formulate wave equations and relativistic heat equations which have the more 
usual properties associated with these equations on a manifold (for example 
they will have a finite speed of propagation). The eigensystem of the edge- 
based Laplacian may be of great use in the study of networks where distance 
and propagation speed are important. The current analysis is currently limited 
to the case of uniform edge lengths. Future work will focus on the case where 
the edge lengths may vary. 
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